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The basis for this research is the 1991 Summer Faculty Fellowship work by the present author. In 
review, Pangia [1992] postulated that the preferred state of a single component infinite plasma is the one that 
will change the least when perturbed. The plasma distribution of such a state would maximize the plasma wave 
damping rate, or, minimizing y where negative y is the damping rate. This explained the tendency for low 
values of k, when fitting plasma distributions to a K-distribution, which is 


(1 + v 2 /(2kv*)) k 

where A K and v K are the overall normalization factor and thermal speed parameter, respectively. A more 
compelling question is why a K-distribution works so well in fitting both electron and ion distributions. 
Hasegawa et ai [1985] derived that the electron distribution subject to a superthermal radiation field is given by 
[1], but that the ion distribution is Maxwellian for speeds less than the electron thermal speed. It is also 
desirable to calculate k, which is not readily doable from the result of Hasegawa et al . [1985]. The present study 
continues from the ideas developed by Pangia [1991] to determine what distribution function maximizes the 
damping rate. A K-distribution will be the outcome for either electrons or ions, with k tending toward 3/2. 

It is necessary to extend the work of the author to include the general case of a multi-component 
plasma. For the sake of completeness, this presentation will reproduce also the arguments of the previous work. 
For a 1 -dimensional, collisionless plasma with electrostatic field, E(t^0, the Vlasov and Maxwell equations are 


afs p. 

9t + V 5x + m s dv ~ 

[2] 

OO 

f = 5>* J dvf s 

[3] 

S -oo 


f =-2>k Jdvvfs 

[4] 


s -oo 


where f s (tpc,v) is the reduced distribution function for particles of type s. In use is the normalization that the 
number of particles of type s that are between x and x+dx with velocities ranging from v to v+dv is f s dx dv. 
Particles of type s have charge qg and mass m s . 

In the absence of any external field, the steady state solution to the Vlasov-Maxwell equations is f s = 
F s (v) and E = 0, where F s is any positive function normalized to the number density. One might then infer that, 
in equilibrium, a plasma distribution is virtually likely to be anything. However, the response of the plasma to 
a perturbation from equilibrium will greatly depend on the function F s (v). Since plasmas evolve to equilibrium 
under the influence of plasma waves, a particular steady state could be favored over all the others based on its 
response to a perturbation. In an actual plasma, waves are usually present, indicating some sort of evolution of 
the distribution function. Although equilibrium may never be fully obtained, the plasma should be at least 
tending toward it. Changes in the background distribution function should diminish as equilibrium is 
approached. Change can possibly be used as the criterion for identifying the preferred steady state of a plasma, 
because a steady state that appreciably changes after having been perturbed is expected to evolve to a different 
steady state which is more stationary. Therefore, it will be postulated that the plasma will tend toward the steady 
state which changes the least when perturbed from equilibrium. 

Since the steady state distribution in a plasma where external forces are absent is spatially uniform, 
each distribution function will be expressed as a sum of a spatially uniform part, fs(t,v), and the deviation from 
uniformity, f ’ s (t,x,v), with a similar expression for the electric field 

f s (t,x,v) = f s (t,v) + f’s(t,x,v) [5] 

E(t,x) = E(t) + E'(t,x) [6] 


XL-1 



The uniform electric field, E(t), may be zero, but, in general, it will exist if a zero wave number mode is part of 
the perturbation or if it develops in time. A spatial average over all space will be used to formally define the 
uniform part of any function Y(x); namely, 




[7] 


The subscript x on the angular bracket is used to identify it as a spatial average. In regard to [5] and [6], f s (t,v) = 
<f s (U,v)> x and E(t) = <E(t,x)> x . Taking the spatial average, [2] through [4] become 


i fs(t ’ v) + { E(t > f ^’ v > + < E ' f ’s>* ) = 0 


*s 

oo 


0 = ^47tqs Jdv f s (t,v) 


S -OO 


-E(t) = - ^4jcqs Jdv v f s (t,v) 


[ 8 ] 

[9] 

[ 10 ] 


Subtracting [8] through [10] from [2] through [4], respectively, gives 

1f + v |r + it* { E(t) f '* + E ' + E ’ f 's - < E ' f ’s>x } = 0 

oo 

3E' 


^=2 4 "qs Jdvf’s 


[ 11 ] 

[ 12 ] 


3F 

dt 


= - ^r4jKk Jdv V f’ s 


[13] 


Equations [8] through [13] are basically the starting equations for quasi-linear theory [see for example, Swanson, 
1989], just extended to include a uniform electric field. Derivable from [8] through [13] are the conservation 
laws for the system, which are written 


oo 

^ Jdv f s (t,v) = 0 

-OO 

oo 

^Tm s Jdv v f s (t,v) = 0 

S -oo 


d 

dt 


I-s 

s 


Jdv v 2 f s (t,v) + rj- E 2 (t) + <E' 2 > x 


\ 

= 0 

) 


[14] 

[15] 

[16] 


Due to the quadratic non-lineanty present in the Vlasov equation, an additional simplifying assumption 
will be made in order to derive a result. Having perturbed a system from equilibrium, it will be assumed that, 
eventually in its evolution, the predominant wave modes that remain are k=0 and k nearly zero, where k is the 
wave number, all other modes having been depleted. This assumption is made even if an instability exists for 
some particular range of modes. The rationale is that non-linear effects will channel wave energy to other modes. 



most of which damp away. To understand the basic effect of the quadratic non-linearity, consider the spatial 
dependence of any mode k at some time L It will have the basic form ± e _ikx . The product of two modes, 
ki and k2, will result in modes ki + k 2 , and ki - k 2 - In particular, if ki = k 2 , there will be mode conversion to 
2kj and 0, and if kj and k 2 are nearly equal, the nearly zero wave mode will result If ki is unstable, the process 
will convert it to stable modes. Modes that find their way to the instability will have been damped beforehand, 
amplified during the instability, and, finally, converted back to stable modes. Therefore, all wave activity will 
experience some damping, either directly or indirectly, during some time in the system’s evolution except for the 
k=0 mode. Due to the initial perturbation having a small amount of total wave energy distributed over all the 
modes, it is being assumed that even an unstable mode would be depleted of its field intensity through the 
eventual conversion of all modes to k=0. There comes a stage in the plasma evolution when only the undamped 
mode, k=0, and the least damped mode, k nearly zero, have survived. 

When the plasma reaches the stage of predominantly k=0 and nearly zero wave modes, further 
amplification of the k=0 mode will be negligible, because only the nearly zero mode, which is decaying, is 
available to convert to k=0. At this point, E(t), which is the k=0 mode, is well-approximated as an oscillation at 
the plasma frequency. From [10], f s (t,v) will oscillate with opposite phase to E(t). Therefore, on a time average, 
the effects of the k=0 mode on changing the distribution function will be zero. Defining a time average for any 
function Y(t) by 


t+T 

<Y> t = jyY(0 [17] 

t 

where T is the plasma period, the time average of [8] is 

| <fs(t,v)> = f ’s>x.t [18] 

where the combined subscript x,t on the angular bracket indicates a double average over space and time. 

The nearly zero mode, which is given by E'(t,x), will be in a regime where its evolution and effects are 
describable by quasi-linear theory. In the quasi-linear approximation, products of E‘ and f' s are retained in [8], 
but neglected in [11]. In search of the distribution that changes the least, a restriction will be made to 
distributions that change only slightly. Such a subclass of distributions can be insured to exist by making the 
perturbation small enough. Specifically, based on [16], the initial wave field must satisfy the condition 


oo 

E 2 (t) + <E' 2 > x « 4n^m s Jdv v 2 f s (t,v) 

s 

at t=0. Then E(t) will be small, and the term E(t) f ' s will be negligible in [1 1], reducing it to 


[19] 


3fs 3fs qs^.dFs 
3t + v 3x + m s dv 


[ 20 ] 


where f s (t,v) was replaced by F s (v) since it only slightly varies. Equations [12] and [20] are the common 
linearized equations for electrostatic waves. The solution is that, for a stable mode of wave number k, E' 
oscillates with frequency to and exponential decays with damping rate hd (negative y corresponds to damping) 
given by [Nicholson, 1983] 


e(©+iyjc) = 1 


47cqj(k 2 m s ) Jdv 
C 


dFs/dv 


= 0 


v - 


[ 21 ] 
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where the integration is in the complex v-plane along the Landau contour, and e is the dielectric function. With 
E' exponentially decaying at a damping of lyl, [18] says that, on the average, the change in f s (t,v) decays at least 
as fast. Consequently, the distribution that changes the least will be the one that damps the nearly zero wave 
mode the most From [21], one seeks the F s that maximizes lyl, or minimizes y. Mathematically, the problem 
at hand is under the classification of Calculus of Variations, where F s is varied slightly and a minimum in y is 
sought. 


The functions about which one varies have to be physical. To insure this, constraints must be 
imposed. The conservation laws, [14] through [16], require that 

oo oo oo 

Jdv F s = n s ; ^Tm s jdv v F s = n P ; m s Jdv v 2 F s = n T [22] 

-oo s — oo 5 -oo 

where, n s , n, P, and T are, respectively, the number density of particle type s, the total number density of all 
particle types, the total plasma momentum, and the total plasma temperature,all of which are constants. The 
existence of these velocity moments imply that the high speed dependence of F s must be less than 1/lvP. This 
puts a limit on damping rate, because, for small k, y is proportional to the dFj/dv evaluated at high velocity 
to/k [Nicholson, 1983]. Therefore, the damping is maximized for the function that has the steepest descending 
slope, which is when the high speed dependence of F s tends to 1/lvP. in regard to [1], k tending to 3/2 would 
maximize damping. 

To find the function for all velocity, additional constraints must be imposed to keep F s physical. One 
condition is that F s must be positive. Another condition arises from the fact that each distribution function F s 
will consist of a definite number of plasma components, A s , which is determined by the origin of the plasma. 
Therefore, it is a fixed property of the system. Defining F s c as the distribution function for component c of 
particle type s, the total distribution for particle type s is 

A s 

F s =£f S iC [23] 

C=1 

The number of components, A s , is determined by counting the humps in F s . By constraining F SiC to be 
positive with only one hump, F s is insured to be physical. This is done by maximizing the following integral 


OO 

Jdv (F s>c ) r 

-OO 


[24] 


where the exponent r has to be determined, where r*l, so that maximizing [24] is independent from the first 
equation in [22]. The problem is now well defined, whereby y is minimized in [21] subject to the constraints in 
[22] while maximizing [24] . A solution for arbitrary and small variations in F SfC only exists if k approaches 
zero. The answer is that each component is given by a K-distribution with appropriate flow speed, where k=1/(1- 
r) with 3/2 < k < «>, and k tending to 3/2 gives the absolute minimum in y . 
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